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ABSTRACT 


A time-independent, linear, plane and axially- symmetric stability 
analysis is performed on a self-gravitating, plane-parallel, isothermal 
layer of non-magnetic , non-rotating gas. The gas layer is immersed 
in a plane-stratified isothermal layer of stars which supply a self- 
consistent gravitational field. Only the gaseous component is perturbed. 
Expressions are derived for the perturbed gas potential and perturbed 
gas density that satisfy both the Poisson and hydrostatic equilibrium 
equations. The equation governing the size of the perturbations in the 
mid-plane is found to be analogous to the one-dimensional time-independent 
Schrodinger equation for a particle bound by a potential well, and 
with similar boundary conditions. The radius of the neutral state is 
computed numerically and compared with the Jeans' and Ledoux radius. 

The inclusion of a rigid stellar component increases the Ledoux radius, 
though only slightly. Isodensity contours of the neutral or 
marginally unstable state are constructed. Large flattened objects 
with masses of 5xl0 6 and radii of 1-2 kpc result. These numbers 
are not inconsistent with the large-scale structure observed in the 
gaseous component of spiral arms in the Galaxy. The possibility is 
discussed that the gravitational instability of the gaseous component 
excites density waves of the type described by Lin. 
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I. INTRODUCTION 

The best-known perturbation analysis with applications in galactic 

astronomy was carried out by Jeans (1928) and appropriately is called 

the Jeans' instability problem. The initial equilibrium state was 

taken to be a self-gravitating infinite uniform gas with a uniform 

gravitational potential field throughout. The last two assumptions 

are, however, inconsistent with the Poisson equation. It is this 

inconsistency which has prompted others to reinvestigate Jeans' analysis 

with self-consistent equilibrium density and potential distributions. 

Jeans applied perturbations to the momentum, continuity, 

Poisson, and heat equation, and after linearization and Fourier analysis 

1/2 

found the condition for marginal stability to be Aj = c s (rr/Gp o ) , 
where c s is the isothermal sound speed of the medium and p Q is its 
density. Any disturbance with length greater than Aj is gravitationally 
unstable and will collapse in a finite time; any disturbance with length 
smaller than Aj is stable. 

The first modification to Jeans' analysis of importance to galactic 

astronomy was undertaken by Ledoux (1951), who considered the more 

realistic initial equilibrium state of a plane-parallel, non-rotating, 

isothermal, self-gravitating gas layer (no stars). The gas density 

at the plane of symmetry assumes its maximum value p Q and decreases 

with distance from the mid-plane according to the well-known formula 
2 

p(z)/p Q = sech (z/H), where H is the scale height of the gas layer 
in the z direction. Using these self-consistent initial conditions, 

1/2 

Ledoux finds the condition for marginal stability to be = c s (2ir/Gp 0 ) , 
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precisely the same relation found by Jeans except that p Q is replaced 
by Pq/2. 

Goldreich and Lynden-Bell (1965) considered the gravitational 

stability of a plane-stratified, self-gravitating, uniformly-rotating 

disk of gas (no stars). They found that (i) pressure effects stabilize 

short wavelength disturbances while rotation stabilizes long wavelength 

2 

disturbances and (ii) when the quantity rrGp/4Q is greater than about 
1.0, where p is the mean gas density and Q is the angular velicity, 
disturbances of dimension several times the layer thickness become 
unstable . 

In what follows, we consider a stability analysis similar to that 
of Ledoux (1951), the principal difference being that the gas layer is 
immersed in a rigid non-perturbable star layer which supplies a self- 
consistent gravitational field. The effect of a magnetic field on the 
stability will be considered in Paper III of this series (Kellman 1972a); 
the effect of a magnetic field plus cosmic-ray gas will be investigated 
in Paper IV (Kellman 1972b). We will not address ourselves directly 
to the grand overall design of spiral structure observed in galaxies, 
and discussed by Lin and Shu (1964, 1966) and Lin, Yuan, and Shu (1969). 
This is probably a collective stellar phenomenon and therefore lies 
outside the range of our analysis. Rather, we limit our attention 
to the existence of the large-scale structure within the gaseous 
component of spiral arms. 

A 21-cm southern survey by McGee and Milton (1964) revealed that 
the principal elements of the gaseous component of spiral arms in the 
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Galaxy are enormous flattened objects, 1-2 kpc in size and 10^ 1^ in mass. 

In addition, these structures were observed to be strung out along the 
length of a spiral arm like beads on a string, their size increasing 
with increasing distance from the galactic center beyond the Sun's 
distance R 0 . Because the projected central density of these structures 
is observed to be only slightly larger than the projected ensity of 
regions between the structures, a linearized stability analysis of the 
type carried out below may shed light on their existence. 

II. STABILITY ANALYSIS 

With these ideas in mind, we proceed to examine the stability of 
a self-gravitating plane-stratified distribution of gas to linear, 
t ime- independent , plane and axially-symmetric perturbations. The plane 
of symmetry is the mid-plane z=0 where the gas and star densities attain 
their maximum values. The axis of symmetry is perpendicular to the 
symmetry plane and arbitrarily positioned. The gas is immersed in the 
gravitational field of a plane-parallel distribution of stars. Both 
components of the two-fluid mixture are assumed to be isothermal, but 
each has a different 'temperature'. We should stress that only the 
gaseous component is perturbed; the stars merely supply a self-consistent 
gravitational field which adds to the gravitational field of the gas. 

This initial equilibrium state is essentially that calculated by 
Kellman (1972c) in Paper I of this series. By restricting the analysis 
to be time-independent, only the marginally unstable or neutral state 
is determined. The inclusion of a stellar component should lend 
the analysis an increased sense of reality as concerns the stability 
of the gaseous component of the Galaxy. 
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The time-independent perturbations may be written as follows: 
p g (r,z) = P e gO) + Ap g (r,z) (1) 

P g (r,z) = P eg (z) + Ap g (r,z) = <v\> P g (r,z) (2) 

cp g (r,z) + cp*(r,z) = cp eg (z) + cp e *(z) + Acp g (r,z). (3) 

p, p, and cp are the density, pressure, and gravitational potential. 

2 

<v tz > is the mean square z turbulent gas velocity dispersion. The sub- 
scripts g and * refer to the gas and star components, respectively. The 
subscript e denotes the equilibrium state described above. A denotes the 
perturbed quantities, functions of both r (the distance from the symmetry 
axis) and z (the distance from the symmetry plane). Since only the 
gaseous component is perturbed, we can write that 


cp*(r,z) = cp e *(z). 


(4) 


so that equation (3) becomes 

cp g (r,z) = <P eg (z) +Acp g (r,z). 


(5) 


Equations (1), (4), and (5) are substituted into the gas hydrostatic 


equilibrium equation 


<v tz > 

— % ' 7(<p g + ^' 


( 6 ) 


and retaining terms only to first order in Ap and Acp we find that 

O O 

O) 


J- V(Ap„) - Z!*S 4p - 1 


eg 


S' 


g 


<vf „> 


8 ' 


eg ' v t;z ' 

Goon (1966) has suggested the usefulness of defining a new 
variable e(r,z): 


e(r,z) = 


Ap g (r,z) 

P e g( z ) 


( 8 ) 
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Equation (7) can now be written in terms of e: 

Ve = — 4— V(Acp ). (9) 

8 

Equation (9) is integrated to give . 

e = + C, (10) 

tz 

o 

where C is a constant. We apply the Laplacian operator V to equation (10) 


with the result that 


V e = <v2 > 
tz 


2, = -V- V 2 (Ac p„) = - ^f^Ap g = _ P eB / PeBO , 


eg' ^ego’ 


where the Poisson equation 


V (Acp g ) = -4rrGAp g 


has been used and where 

8 \ 


is a scale height for the gas in the absence of stars. Equation (11) 
may be solved by the method of separation of variables, in which case 
e is written as 

e (r , z) = X(z) Y(r) . (14) 


After some simplification, equation (11) becomes 


_ (Y"(r) + 7 Y*(r)) _ 2 _ (X"(z)+X(z) l/2S? (z)/p pf J 

K. 


g e g 


X(z) 


ego 7 (15) 


The separation constant k arises because the left-hand side of 
equation (15) is a function only of r while the right-hand side is a 
function only of z. 

The differential equation for Y(r) 


Y"(r) + 7 Y’ (r) + k 2 Y(r) = 0 


(16) 
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has the solution 

Y(r) = J Q (kr), (17) 

where J Q is the zero order Bessel function. The other linearly indepen- 
dent solution to equation (16) -* ■» as r -* 0 and is therefore not 
considered further. The differential equation for X(z) = X k (z) is 

X k (z) +( 2Hg P eg (z)/p ego" k2 ) x k< z > = °» (18) 

with the boundary conditions 

X£(0) =0 (19) 

lim X k (z) = 0. (20) 

z |-* 00 

In addition, X k (0) is normalized to unity. Equation (19) arises 

because Apg(r,z) is an even function of z. Equation (20) expresses the 
constraint that Ap g (r ,z)/p e g(z) — 0 as |z| -* 00 . 

The general solution to equation (11) can be expressed as an 
integral over k of the product X^(z)J Q (kr) with the appropriate expansion 

amplitudes A(k) 

e(r,z) = A(k)X k (z)J 0 (kr)dk, (21) 

o 

where J Q (kr) is the solution to equation (16) and X k (z) is the solution 

to equation (18) subject to the boundary conditions imposed by equations 
(19) and (20). Equations (10) and (21) can be combined to give an 
expression for the perturbed gas potential: 

Acp g (r,z) = <v£ z > | A(k)X k (z)J 0 (kr)dk - <v\^> C. 


( 22 ) 
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The solution of the Poisson equation for the perturbed gas potential 
is just the sum of the particular solution of the inhomogeneous equation 
and the general solution of the homogeneous equation. Equation (22) 
with C=0 corresponds to the former. The general solution of the homo- 
geneous equation is equivalent to the general solution of the Laplace 
equation 

V 2 Ap gh = 0 3 (23) 

where the subscript h denotes homogeneous. 

To solve equation (23) we expand Acp gk in terms of J 0 (kr) with 

expansion amplitudes s^(z): 

Acpg h ( r 3 z) = <v^ 2 > s k (z) J Q (kr)dk. (24) 


Expressed in terms of cylindrical coordinates, equation (23) assumes 
the form 


1 i_ 

r Sr 


(r fj A<p gh ) + 



Acpgh = 


0 . 


(25) 


Equation (24) is substituted into equation (25), and after some 
differentiation and algebraic manipulation we find that 

,(kr) 


i 


s k (z) J Q (kr) + ks R (z) 


( 3j r 


+ kJ (kr) 
o ' 


dk = 0. 


(26) 


We can make use of the well-known relations between the zeroth and 
first order Bessel functions J Q and J-^ and their derivatives 


J Q (kr) = -J 1 (kr) 

Jg(kr) = -J Q (kr) + ^ J!(kr) 


(27) 

(28) 


to simplify equation (26): 
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0 = (s^z) J Q (kr) - k 2 s k (z) J Q (kr))dk. (29) 

Since the J Q functions are orthogonal, it follows that 

sj^(z) - k 2 s k (z) = 0, (30) 


which has the solution 
-k| z 


s k (z) = w 


lk e 


+ w 


+k z 


2k fe 


(31) 


The coefficients W 2 k must vanish if the boundary condition lim Ac p 

UI-HCO 


gh 


= 0 


is to be satisfied. 

We have now obtained the functional form of s k (z) in terms of which 
Acf>g k has been expanded. Substitution of equation (31) into equation (24) 
easily gives 

Acpg h (r,z) = <v\^> J q w lk e k l Z l J Q (kr)dk. (32) 

Note that Acpg k (r,z) expressed here satisfies the boundary condition that 


the perturbed gas potential be symmetric with respect to the symmetry 
plane z = 0. Adding equations (22) (with C = 0) and (32) then gives the 
general solution to the Poisson equation satisfying, of course, the gas 
hydrostatic equilibrium equation: 

Acp g (r,z) = <vj; z > (A(k)X k (z) + w^e k ^) J Q (kr)dk. (33) 

For the sake of completeness we may easily derive from equations 
(8) and (21) an expression for the perturbed gas density satisfying 
the Poisson and gas hydrostatic equilibrium conditions: 

APg( r ,z) = P eg (z) f Q A(k) X k (z) J Q (kr)dk. 


( 34 ) 
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It is appropriate here to consider more closely the nature of 
equations (18), (19), and (20) as regards the nature of the marginally 
unstable state. Equations (18) - (20) define what is known as the Sturm- 
Liouville problem from the theory of differential equations. This is 
the problem of determining (a) the relation between the separation constant 
k and the function X k (z) and (b) the influence on k of the boundary conditions 

imposed on X k (z). Stated more simply, we are searching for those values 

of k that lead to X k (z) satisfying both equation (18) and the boundary 
conditions imposed by equations (19) and (20). It is worthwhile to note 
the mathematical similarity of our problem and the quantum mechanical 
problem of determining the allowed energies E of a particle bound by a 
one-dimensional potential well V(z). The particle wave function 'F(z), 
analogous to X k (z), satisfies the one-dimensional time-independent 

Schrodinger equation similar to equation (18), and is subject to the 
square infcegrability condition that | Y 2 ( z ) | d z be finite, similar to 

2 

equation (20). To carry the comparison one step further, k corresponds 

2 

to the particle energy E and 1/2H p (z)/p corresponds to the form 

g eg ego 

2 

of the potential well V(z). Only certain discrete values of k will 

yield X^Cz) that obey the boundary conditions. Since Ap (r,z)/p (z) = 

8 ® 8 
co 

J q A(k) X k (z) J o (kr)dk, the statement X R (z) - 0 as |zj is equivalent 

to the statement Ap g (r ,z)/p eg (z) - 0 as |z| - oo. it is the latter 

condition then that makes this an eigenvalue problem. 

The stellar component will make its presence felt only as it 

2 

effects the depth and or width of the potential well 1/2H o (z)/o 

g K eg v ' 'ego’ 


k 
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or equivalently, as it effects the distribution of gas p (z) above the 

eg' 

galactic plane. It is at this point then that the present stability 

analysis differs from the earlier work of Ledoux (1951). Ledoux considered 

an isothermal layer of gas (no stars), so that the density distribution 

2 

obeys the well-known relation p eg (z)/p eg0 = sech ( Z / H )* When a stellar 

component is included, P eg (z)/p ego can be found by solving simultaneously 

the gas and stellar hydrostatic equilibrium equations and the Poisson 

2 1/2 

equation and by choosing appropriate values for p ego , p^ 0 , <v£ z > , and 

2 1/2 

<v^ z > . This has been done in Paper I, from which we recall equation (11): 


in p (z)/p 

dz z M eg' K ego 


4ttG 




p (z) + p (p (z)/p ) <v tz > / <v *z > 

eg v ' K *o vp eg v " ego 7 


solution to equation (18) in the limit z -» + is just 

(+) -kz (+) +kz 
X k (z) = <* k e + 3k e • 


(36) 


(35) 


The observed values p egQ = 1 H atom/ cm 3 = 0.025 %/pc 3 (Weaver 1970), 

- 2 1/2 

p^ 0 = 0.064 P^/pc (Luyten 1968), and <v^ z > = 18 km/sec (Woolley 1958) 

2 1/2 

are chosen, <v tz > is allowed to vary between 0 and 20 km/sec. 

At this point we digress slightly to discuss the numerical approach 
employed in computing the discrete value or values of the separation 
constant k. It can be shown (see Appendix, equation (A16)) that the 


We wish to find the restriction on k and X k (z) if P^^is to vanish. 


If Pi 


(+) 


= 0, it is clear that 


X k (z) + kX k (z) = 0. 


(37) 


To derive the eigenvalue k from the restriction dictated by equation (37), 
a sufficiently large z is considered, denoted by z Q , 


for which equation 
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(36) holds with good accuracy. Employing numerical methods, equations 
(18) and (35) are solved simultaneously for X^(z), allowing us to tabulate 

' i 

Xk(z c ) and X^(z c ) • The quantity P = Xk(z c ) + kXk(z c ) is then calculated 

at each of four equally spaced points covering the range of k. The 
eigenvalue k lies in the interval over which P changes sign. This 
interval is divided into three subintervals and the process repeated 
until the desired accuracy is attained. 

To determine how the eigenvalue k relates to the linear size of the 
neutral state , we recall equation (17) which expresses the solution to 
the second order differential equation for Y(r): 

Y(r)=J 0 (kr). (17) 

We note that k must have units of reciprocal length, and that for a 
given value of k the first zero in Y(r) occurs when 

r l = a l^ k l’ ( 38 > 

where is the first zero of the zero order Bessel function J Q and k-^ 

is the smallest eigenvalue found by the numerical methods enumerated 
above. More generally, if more than one eigenvalue k is found, we 
have that 

r n = °n /k n’ < 39 > 

where o? n is the nth zero of J Q . Focusing on equation (38) we observe 

that k-^ determines the radius of the 'perturbation', or more accurately 

the radius r^ of the marginally unstable state in the plane of symmetry 


z=0, since 
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e(r! = o^/ki, z=0) = X k (o) J Q (<*]_) = A P g (^ 1 / k 1 , 0)/p eg ( 0) = 0 . (40) 

When equations (18) and (35) are solved simultaneously subject to 

the boundary conditions imposed by equations (19) and (20) and subject 

to the values chosen for p , p , <v^ and <v^ only 

ego’ *o *z ’ tz 5 J 


one 


2 1/2 

eigenvalue k is found. Since <v tz > is allowed to vary between 0 and 

2 1/2 

20 km/sec, k. is determined as a function of <v_ > . r, may be found as 

•L tz 1 J 

2 1/2 

a function of <v tz > from the relation rj_ = a^/ki. The results are 

2 1/2 

plotted in Figure 1. We see that r^ is directly proportional to <v tz > 

2 1/2 

and increases from 1 to 2 kpc as <v tz > varies between about 8 and 
16 km/ sec . 

2 1/2 2 1/2 

It is useful to inquire how p , p , <v > , and <v > effect 

ego’ *o’ *z ’ tz 

kj_ and thus r-^. Recall that the equivalent potential well has the form 

2 ( 2 
V(z) = 1/ 2Hg p eg (z)/p ego . ' The central depth is just V(0) = l/2H g since 

2 

P e ~(0)/ p. n = 1. Since H is a function only of p and <v >, p and 

'"6 o G§0 t Z /VO 

2 

<v > have no effect on the central depth. Increasing the depth of the 

2 

potential well by decreasing H 2 (and thus increasing p or decreasing 


2 

<v tz > ) will increase k^ and decrease r-^. p and <v z > will, however, 

effect the width of the potential well, in the sense that increasing 
P or decreasing <v > decreases the width which causes ki to decrease „ 

VrZ J- 

and r-^ to increase. In summary, the depth is determined only by p e g 0 
2 

and < v tz >, while the width is determined primarily by p , and to a 

2 2 

lesser extent by P e g Q , <v* z >, and <v tz >. The inclusion of a stellar 
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component will increase the radius r^ of the marginally unstable state 


above its value when p = 0. 

*o 


It is instructive to compare r-^ calculated with P e g Q = 0.025 1^/pc , 


o 2 1/2 

P*o = 0.064 1^/pc , and <v^ z > = 18 km/sec with the Jeans' radius 

(\ t / 2) and the Ledoux radius (X T / 2 ) , each calculated with p = 0.025 
J L ego 

/ 3 2 1/2 . / 

%/pc . < ' v t z > I s a free parameter, allowed to vary between 1 and 20 

km/sec. The results are presented in Table 1. It is particularly useful 

to compare r-^ (Kellman) to r-^ (Ledoux), since both assume plane-parallel 

gas distributions in the equilibrium state. The presence of a rigid 
stellar component with p^ D = 0.064 lY^/pc 3 and < v 2 z > 1 ^ 2 = 18 km/sec 

increases the radius of the marginally unstable state, though only slightly. 

Isodensity contours of the marginally unstable or neutral state 
may be derived by combining equations (1) and (34), with the result that 


Pg(r ,z) ^ P eg (z) 


(1 + A(k) X k (z) J Q (kr)). 


(41) 


ego 


ego 


The integral over k in equation (34) reduces to just one term since 
only the neutral state is being considered and since only one eigenvalue 
was found. At the center of the perturbation r = z = 0, X^(z) = 1, 

J Q (kr) = 1, and p eg (z)/p ego = and equation (41) becomes 


p g(°> 0 ) . , ... N 

— S = 1 + A(k). 

p ego 


(42) 


21-cm observations of large gas structures in spiral arms (McGee and 
Milton 1964) indicate that 1.0 S p (0,0)/p 

O 


^ 2.0, and therefore 
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0 £ A(k) fi 1.0. Making use of equation (35) for p (z) / p , equation 

eg ego 

(18) for X,(z), and numerical tables for J Q (kr), equation (41) is used 

K 

to construct an (r,z) grid of values for p (r,z)/p from which 

g ego 

contours of equal p g (r,z)/p eg0 may be derived. Figure 2 displays three 
isodensity contours p g (r,z)/p eg0 = 1.3, 0.7, and 0.3 calculated with 

A(k) = 0.5. The neutral state is clearly quite flattened, and is 
reminiscent of contours presented by McGee and Milton (1964) and 
Westerhout (1956), the latter referring to cross sections of spiral arms. 

The mass of the marginally unstable state may be expressed as an 
integral of the gas density over a cylinder perpendicular to the galactic 
plane, of infinite height and of radius r^ = a^/ k^: 


M g = 2n J ^ p g (r , z)r dr dz. 


(43) 


Equation (43) can be evaluated numerically with the aid of equation (41). 

2 1/2 

Choosing <v tz > = 10.0 km/sec (see Paper I) and values quoted previously 

2 1/2 6 
for p , p , and <v > , we find that M„ » 5x10 M^. 

'ego’ H *o’ *z g ^ 

Equations (18) and (35) may be written in dimensionless form 

^2 % - - < | - c k 2 > x k < 44 > 


Ivot = - - (a + ya s ) , 
dp 2 2 


(45) 


where the dimensionless quantities C^, a , 0 , 7 , and 6 are defined 


as follows: 
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c k = kH g 

(46) 

a ~ Peg^ z ^Pego 

(47) 

3 = 

(48) 

^ P *o /,p ego 

(49) 

6 = <\>/<\>- 

(50) 


The radius of the neutral state in the plane of symmetry can be 
written as r-. = but since H = (<v^ >/8nGp ) 1 ' ,2 ,r. depends 

2 1/2 

on the ratio (<v >/p ) , precisely the functional dependence found 

v tz ego' 

2 1/2 

by Jeans and Ledoux. If <v tz > is independent of distance from the 


galactic center beyond the solar distance R Q (see Paper I), r-^ would be 


expected to increase with increasing distance from the galactic center 
since P e g Q decreases with increasing R. Westerhout (1956) finds that 

(corrected to the new galactic distance scale) P ego ( R = H k P c )/P ego ( R = 15 kpc) 


ss4.7 from which it follows that r-^(R - 15 kpc)/r-^(R - 11 kpc) 

ss 2.2 due to a change in p alone. In conclusion, since ri = 2 _ -4 _ 0 5 ji 

ego c k S 

2 1/2 

r . depends principally on <v > and p through H to a smaller 
1 r r tz ego g 

2 1/2 

extent on p through C, , and to an even smaller extent on <v_ > 

tc z 


through C k . 
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III. COMPARISON WITH THE OBSERVATIONS 
It is interesting that a 21-cm southern survey by McGee and Milton 
(1964) revealed that the principal elements of the gaseous component of 
spiral arms in our galaxy are large flattened structures with dimensions 
10 7 and 1-2 kpc , strung out along the length of the arms much like 

beads on a string. More recently, Kerr (1968) and Westerhout (1971) 
have pointed to the existence of large gas condensations forming the 
major components of spiral arms. The size and mass of the individual 
gas structures observed by McGee and Milton are close to the values 
we have calculated for the marginally unstable state. Further, their 
linear size is observed to increase markedly with increasing distance 
from the galactic center beyond R Q , consistent with our findings above , 

while their mass is observed to remain essentially constant at about 
10 7 N^t). In addition, they occur predominantly in the outer arms beyond 
R Q (McGee 1964). It is also of interest that McGee and Milton (1966) 

have observed similar gas structures in the Large Magellanic Cloud with 
a mean linear size and mass of about 600 pc and 4x10 respectively. 

2 1/2 

The smaller size may indicate that P e g D is larger, <v tz > is smaller, 
and or p is smaller than corresponding values in our galaxy, if the 
theory is to be believed. 

IV. DISCUSSION 

Recent work by Lin (1970) has refocused attention on the importance 
of classic Jeans' type stability analyses of the gaseous component of 
galaxies. It is Lin's opinion that density waves (primarily in the 
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stellar component) are initiated or excited by the gravitational 
instability of the gas in the outer regions of galaxies, where the 
star density no longer overwhelmingly dominates the gas density and where 
the gas turbulence, which tends to inhibit instability, may be small 
because of reduced energy inputs. Excited in the outer regions, the 
density waves travel inward, subsequently giving rise to two-armed 
spiral patterns. Once formed, the large gas condensations are stretched 
by differential galactic rotation into trailing structures, and because 
they lie near or outside the corotation distance, form material arms. 

The pattern speed Qp over the whole galactic disk is determined then, if 
Lin is correct, by the angular velocity of these outer gas condensations 
around the galactic center. In support of Lin's ideas, Shu, Stachnick, 
and Yost (1971) obtained a good fit to the spiral structure observed 
in M33, M51, and M81 by assuming that Qp is equal to the angular rotation 
speed of the outermost HII regions. 

We mentioned above that the large gas condensations observed by 
McGee and Milton (1964) are preferentially located at distances beyond 
the solar distance R_ . Not only does this fact agree nicely with Lin's 

ideas, but it follows directly from results obtained in Paper I of this 

series. There we found evidence that the rms z turbulent gas velocity 

2 1/2 2 2 1/2 
dispersion <v tz > (or more correctly Q = (<v tz > + B Q /8rrp eg0 + P c - r0 ./P ego ) ) 

increases with decreasing distance from the galactic center in the range 
4 kpc ^ R s 10 kpc . It follows then that gravitational instability of 
the gas disk is inhibited when R < R Q . 
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APPENDIX 

ALTERNATE CALCULATION OF Ap (r,z) 

O 

As a check on the correctness of equation (22), Acp (r,z) may be 

§ 

derived by a different approach, beginning directly with the solution 
to the Poisson equation for the perturbed gas potential. Expressed in 
cylindrical coordinates, we have that 

«„■ 1”1 'V 1 ' "■ z > rdc <“ Wz - < A1 > 

(r Q , 0 Q , z Q ) and (r, 0 , z) are the coordinates of any two points and 
R is the distance between them. It can be shown that 

R = m S 0 e m cos t™( e - 0 ort J 0 J m< kr > J m< kr o> e ^ Z °' dk > ( A2 ) 

and noting that 
„2tt 

J cos[m(0-0 o )]d0 = 2n6 mo , (A3) 


-k z-z 


ol dk. 


where 6 is the Kronecker delta symbol, it follows that 

Acp g (r o’ z o ) = 2x10 So rdr J-o= dz A Pg< r ’ z > Jo V kr)J o (kr o> e 

The perturbed gas density APg(r,z) may be expanded in terms of the zero 
order Bessel function J 0 (k'r) and the appropriate expansion amplitudes 
T k ' (z) 

CO 

A p g (r,z) = J o T k ,(z) J Q (k'r) dk' , (A5) 

and if we note that 

J Q r V kr > Jo( k ' r ) dr = { S(k'-k), (A6 ) 

it follows that 

CQ °C I. 

Acp g (r 0 ,z 0 ) = 2rrG J^dz T k (z)e" k| Z ' Z ° 1 J Q (kr)dk. (A7) 


(A4) 
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Note that two expressions have been derived for Ap g (r,z), represented 
by equations (34) and (A5) . If they are equated we find that 

T k (z) = A(k)p eg (z)X k (z) , (A8) 

and equation (A7) becomes 

Ac P g ( r o, z 0 ) = 2 ttG J o dk . A ,(_ k .)J.o( kr o) p eg (z)X k (z)e" k ^" z o!dz. (A9) 


Equation (18) can be recalled to supply an expression for the quantity 
p eg (z)X k (z) appearing in equation (A9): 

Peg(z)X k (z) = 2p ego Hg (k 2 X k (z)-X k (z)). (A10) 

It is this condition that, when substituted into equation (A9), constrains 
Acpg(r Q , z Q ) to satisfy the gas hydrostatic equilibrium equation: 


Acp g (r 0 , z Q ) = 4 ttGH, 


'g H eg 


cj'l „ k A(k)J 0 (tro)j ^ ( k 2 X k( z ) -X k( z))e‘ k ^‘ Z °ldz. 


To proceed further it becomes necessary to evaluate the integral 
over z in equation (All). This can be accomplished by splitting the 
range into two parts: 


J'l (k 2 X^(z)-X^(z))e" k Z_Z ° dz - r f° 


+ k 2 J*j X k ( z)e +k < z - z »>dz - j; o>a, + k 2 f 


The first and third terms on the right hand side of equation (A12) are 
integrated twice by parts and added to the second and fourth terms, and 
after some simplification we find that 

' /„\ + k ( z " z n> 


J (k 2 X k (z)-X k (z))e ° dz = -X k (z)e 


o 

-oo 


+ kX k (z)e 


+k(z-z 0 ) 


t s - k(z-z 0 ) 


kX k (z)e 


-k(z-z Q ) 


(A13) 


(All) 


(A12) 
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To evaluate these terms in the limit z - + » we recall the differential 
equation for X k (z): 

V z) + (^-" Pe § (z)/p ego - k2 ) X k< z > = °- (A14) 

In the limit z -» + 00 , equation (A14) becomes 

xj^(z) = k 2 X k (z) , (A15) 

since Mm _1 — P e g( z )/pe 20 = Tke solution t0 equation (A15) 

z * i -00 ® 

~ g 

can easily be written: 


<+> e' k2 + 6k (+> e +kZ 

(A16) 

( ' ) e- kz + B k ( ' ) e +kz . 

(A17) 


' k 

Referring back to. equation (34) for the perturbed gas density 

A P g(r ,z) /p e g(z) = J*" A(k)X k (z)J o (kr)dk, (A18) 

it is clear that unless Mm X, (z) = 0, the quantity Apg(r ,z) /p e g(z) will not 

Z“* + oo k 

(-) (+) 

-» 0 as z -» + 00 . This requires that c* k = P k =0. Symmetry of 

Ap (r,z)/p e „(z) about the plane z = 0 results in the further restriction 
that a,^= Equations (A16) and (A17) therefore become 


z > 0: X, (z) = ye 


-kz 


z < 0: X k (z) = ye 
where y = ar k ^ = 3 k ^ ^ 


+kz 


(A19) 

(A20) 


. It follows that 
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z 


z 


I “ 1c 

> 0: X k (z) = ~ 7 ke 

i -\-\tz 

< 0: X k (z) = /ke 


(A21) 

(A22) 


We can now evaluate each term in equation (A13) at the infinite limit. 
For example, consider the first term: 

lim (-x' (z)e +k(z " Zo) ) = Aim (- 7 ke‘ kz °e +2kz )= 0. (A23) 

Similarly, the other three terms evaluated at the infinite limit also 
vanish. We are now able to evaluate the integral over z appearing in 
equation (All): 


- X k(z))e _k Z ” Z ° dz = -X k (z 0 )+kX k (z 0 )+X^(z 0 )+kX k (z 0 ) = 2kX k (z Q ) 


(A24) 


Equation (All) therefore simplifies to 

Acp g2 (r,z) = SrrGHg p figo A(k)X k (z)J 0 (kr)dk, (A25) 

where the arbitrary coordinates (r Q , z 0 ) have been replaced by (r,z). 

This is, however, just the particular solution to the Poisson equation 
(and consistent with the gas hydrostatic equilibrium equation) , for which 
we have derived another expression in the way of equation (22) 

(with C = 0) : 

Ac P g i( r »z) = <v 2 > A(k) X k (z)J Q (kr)dk. (A26) 

2 2 

Recalling that = <v tz >/8rTP e g 0 , it is clear that 

Ac P g l(r,z) = Acp g2 (r,z). (A27) 

The purpose of the above derivation was to obtain an independent 
check on Acp g ^(r,z) derived in Section II, and the equality expressed by 

equation (A27) gives us confidence that the expression obtained for 


Acp (r,z) is correct. 
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TABLE 1 

COMPARISON BETWEEN THE JEANS', LEDOUX, AND 

2 1/2 

KELLMAN RADIUS AS FUNCTIONS OF <v„ > 

tz 


. 2 1/2 
<v tz > 
(km/sec) 

r-^ (Jeans) 
(kpc) 

r y (Ledoux) 
(kpc) 

r^(Kellma: 

(kpc) 

1.0 

0.085 

0.120 

0.127 

2.5 

0.213 

0.301 

0.318 

5.0 

0.425 

0.601 

0.635 

7.5 

0.638 

0.902 

0.953 

10.0 

0.850 

1.202 

1.270 

15.0 

1.275 

1.803 

1.905 

20.0 

1.700 

2.404 

2.540 
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FIGURE CAPTIONS 

1. The dimension of the marginally unstable state in the symmetry 
plane (z = 0) as a function of the rms z turbulent gas velocity 
dispersion. 

2. Isodensity contours pg(r,z)/p e g Q = 1.3, 0.7, and 0.3 of the 
marginally unstable state, calculated with A(k) = 0.5. 
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